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The SLUGGS survey: Probing the supermassive black hole 
connection with bulges and haloes using red and blue globular 
cluster systems 
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ABSTRACT 

Understanding whether the bulge or the halo provides the primary link to the growth of super- 
massive black holes has strong implications for galaxy evolution and supermassive black hole 
formation itself. In this paper, we approach this issue by investigating extragalactic globular 
cluster (GC) systems, which can be used to probe the physics of both the bulge and the halo 
of the host galaxy. We study the relation between the supermassive black hole masses (Mbh) 
and the globular cluster system velocity dispersions (ctgc ) using an updated and improved 
sample of 21 galaxies. We exploit the dichotomy of globular cluster system colours, to test 
if the blue and red globular clusters correlate differently with black hole mass. This may be 
expected if they trace the potentially different formation history of the halo and of the bulge 
of the host galaxy respectively. We find that A/bh correlates with the total GC system velocity 
dispersion, although not as strongly as claimed by recent work of Sadoun & Colin. We also 
examine the A/bh — ^gc relation for barred/bar-less and core/non-core galaxies, finding no 
significant difference, and for the first time we quantify the impact of radial gradients in the 
GC system velocity dispersion profile on the A/bh — cqc relation. We additionally predict 
Mbh in 13 galaxies, including dwarf elliptical galaxies and the cD galaxy NGC 3311. We 
conclude that our current results cannot discriminate between the bulge/halo scenario. Al- 
though there is a hint that the red GC velocity dispersion might correlate better with 1/bh 
than the blue GC velocity dispersion, the number statistics are still too low to be certain. 
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1 INTRODUCTION 

Extragalactic globular clusters (GCs) may provide key insight into 
the connection between galaxies and superma ssive black holes 
(SMB Hs). GCs are typically old (> 10 Gyr, ISrodie & Straderj 
120061) and may have witnessed the events which formed the 
SMBH in the first place. Moreover, GC systems usually come in 
two subpopulations, t hought to be the result of different forma- 
tion mechanism s (e.g. I Ashman & Z epf 1992; Forbes et al. 1 9971 : 
ICote et al.lll998l) . The blue ( metal-poor) subpopulation has been 
associated with ga l axy halos i Forte et al.ll2005c iMoore et al .11 20061 : 
iForbes et al.ll2012l : ISpitler et al.ll2012h . It may have originated in 
metal-poor dwarf galaxies at high redshift consequently accreted 
into the halo of larger systems ( jElmegreen et al.i,2012i) . The prop- 



erties of the red (metal-rich) GCs are similar t o those of the galaxy 
bulge JStrader et al.ll201 ll : iForbes et al.ll20I2h , perhaps because of 
a coeval format ion, such as in a turbul ent disk dShapiro et al.l20ld) 
or in a merger JKruiissen et al.ll20I2h . Therefore, if the growth of 
SMBHs is primarily driven by recent merger events, one might ex- 
pect a stronger correlation between red GCs and SMBHs . Con- 
versely, if the properties of SMBHs were set during the primor- 
dial formation of their hos t galaxies, we might expect a stronger 
correlation with blue GCs JOmukai et al.ll2008l : JMaver et al.ll2010l : 
lDebattistaetal.ll2013h . 

There exists a surprisingly good correlation between the to- 
tal number of GCs (both blue and red) per galaxy (Nqc) and 
the black hole mass of galaxies (A/bh). However, this does 
not necessarily imply a primary correlation between GCs and 



© 2012 RAS 



Pota et al. 



SMBHs l |jahnke&Maccidl201lh . In fact. I Snyder eTal] 1201 ih ar- 
gued this correlation to be indirect as expected if i t was a conse- 
quence of the debate d black hole fundamental plane JHopkins et alj 
I2OO7I : lGrahanill2008l) . Nevertheless, the Mbh - A^gc relation has 
been shown to have an in t riguingly small scatte r at fixed Mbh 
( iBurkert & Tremainell2oTol iHarris & HairislbOllb . iRhodel (|2012|) 
has recently shown that these findings are driven by low number 
statistics, and that an improved galaxy sample returns a scatter at 
fixed Mbh which is larger than previously inferred. Rhode addi- 
tionally found similar slopes and scatters for the relations for the 
blue and the re d GCs. 

Recently. ISadoun & ColinI 1 120 121) (hereafter SC12), have ex- 
amined the correlation between the GC system velocity dispersion 
and Mbh for twelve galaxies, including the Milky Way. Their re- 
sults suggest a tight correlation between A/bh and the velocity dis- 
persion for both the red and blue GC subpopulations, with an in- 
trinsic scatter e always ^ 0.33 dex, indicating a very tight corre- 
lation. They also find that the red GCs are more closely correlated 
(e = 0.22 dex) with Mbh than the blue GCs (e = 0.33 dex). 

In this paper we revisit the work of SC12 with an expanded 
sample of 21 galaxies and updated A/bh values. We supplemented 
our sample with high velocity res olution data from the ongoing 
SLUGGS survey jPotaet al.ll2013h and we re-analysed literature 
data with the same method. We tested if the tight correlation seen 
for the red GCs is real or driven by sample selection or method- 
ology biases. The outline of the paper is as follows. We describe 
the data in Section |2] and their analysis in Section |3] Results are 
then presented and discussed in Section |4] Conclusions are given 
in Section[5] 
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Figure 1. Difference between the rotation-subtracted velocity dispersion 
a"GC and the rotation-included velocity dispersion vqq without any colour 
split. Filled points and empty squai'es ai'e the data from the SLUGGS suivey 
and the literature respectively. The two quantities are generally consistent 
with each other, but they disagree by up to ~ 40 km s~^ in systems with 
significant rotation. 



of which (NGC 4486) was already discussed in SC12. We use 
the most recent black ho le m ass measurements as sum marized in 
iMcConnell & M j ( l2012h and ICraham & ScottI ( I20T2I) . This gives 
us a sample of 21 galaxies, nearly double the number used by SC 12, 
which are listed in Table I All 



2 GALAXY SAMPLE 

We study a subset of galaxies with direct A/bh measurements and 
with more than ten GC radial velocity measurements. From the lit- 
erature, we compiled a list of 13 galaxies. This includes all the 
galaxies discussed in SCI 2, e xcluding the Milk y Way, and two 
additional galaxie s: NGC 253 JOlsen et alj[2004l) and NGC 3585 
jPuzia et al.ll2004h . not studied by SC12 because the uncertainties 
on the GC velocity dispersion were not quoted in the parent papers. 
T he Mi l ky Wa y is not included in this study because the results 
of lCotg i ll 9991) suggest that the still uncertain velocity dispersion 
of the Milky Way GC system is unusually large for its black hole 
mass. Moreover, the fact that the Milky Way GC analysis is carried 
out in three-dimensions rather than in projection, makes the com- 
parison with other galaxies not straightforward. We also update the 
GC catalo gue used by SC12 for N GC 4594 with the latest com- 
pilation of lAlves-Brito et all 1I2OI ih . We note that SC12 used A/bh 
values from lGiiltekin et al.l ( 120091) although more recent A/bh were 
sometimes available. 

For NGC 224 ( M31) we use the GC system velocity dispersion 
measurements from lLee et al. (2008), because their catalogue is not 
available on-line. 

In reg ard to NGC 253, there are tw o p ublic GC catalogue s for 
this galaxy: iBeaslev & Sharpie^ ( I2OO0I) and lOlsen et aLl ( l2004l) . for 
a total 38 GCs. However, we were unable to find a reliable calibra- 
tion offset between the radial velocities of the four GCs in common 
between these two datasets. We decided to use the Olsen catalogue 
only, because it is larger in size (24 GCs) than Beasley's dataset 
(14 GCs). 

The biggest strength of our data set is the addition of a fur- 
ther 9 new early-type galaxies from the SLUGGS survey, one 



3 METHOD 

3.1 The globular cluster system velocity dispersion 

The stellar velocity dispersion, (Jt , used in the A/bh — cr* relation is 
usually defined either as the luminosity-weighted velocity disper- 
sion within l/S*'' of an effective radius Re, or within 1 Re (ere), 
and/or as the central velocity dispersion (ao). Although they rep- 
resent physically distinct quantities, Oe and gp have b een reported 
to be consistent with each other CGiiltekin et al.ll2009l) . This stems 
from the fact that the velo city dispersion profile s vary only weakly 
within these regions (e.g.. lEmsellem et al.ll201lh . 

The detection of extragalactic globular clusters occurs pre- 
dominantly at /? ^ Re. Therefore none of the stellar velocity dis- 
persion quantities are directly recovered with GC data. We define 
the GC system velocity dispersion in two different ways, which are 
similar to the quantities used for stellar data. This also takes into ac- 
count that some GC systems can have a r otation component w hich 
is as large as that of the random motions jBeaslev et al.l200q) . 

The first quantity, ctgc, assumes a Gaussian velocity distribu- 
tion and it is de fined as the stand ard deviation with respect to the 
model function dCote et al.ll200ll) : 



v{e) 



-I- Urot sin(6'o - 0), 



(1) 



which measures the GC rotation amplitude iirot as a function of 
the azimuth 6, with So being the direction of the angular mo- 
mentum vector and Vsys being the systemic velocity of the host 
galaxy. We use a varia tion on Equation [T] originally designed by 
iKrainovic et alj ( 120061) f or IFU data-cubes a nd then extended to 
sparsely sample d data by Proctor et al.l(l2009h . We then minimise a 
X^ function (see lBergond et al.l2006l) to compute the best fit param- 
eters (wrot, CGc, do). Uncertainties were derived by bootstrapping 
the sample 1000 times to derive 68 per cent confidence intervals. 
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Figure 2. Black hole mass as a function of rotation-included GC system velocity dispersion. Left, central and right panels show the A^bh — ^GC relation for 
all, red and blue GCs respectively. Data from the literature and from the SLUGGS survey are shown as open squares and filled points respectively. The black 
solid line is the best fit to the A'/bh — '^GC relation. The dashed lines are the best-fit to the A/bh — ""GC relations (whose datapoints are not plotted here for 
clarity). The slope and the intercept of the best-fit lines are the average between the values from the forward and inv erse regression (see T ableQ]. The dotted- 
orange line is the stellar M-qh — a, relation from the average between the forwai'd and the inverse regression from lGraham et al.l ( 1201 lb : a = 8.14 it 0.05 
and^ = 5.54 ±0.40. 
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Table 1. Solutions to log(AfBH/A/0) = a + /31og(tTGc /200 km s^^) for different GC subsets. Shown are the sample size A*^, the intercept a, the slope /3, 
the intrinsic scatter e and the total rms scatter A in the log A'/bh direction for both the forward (minimise log A/bh residual) and for the inverse regression 
(minimise log a residual). 



We will refer to the rotation-subtracted velocity dispersion of the 
red, blue and all GCs as ctgc.r, ctgcb, ctgc respectively. 

The second quantity, ugc. does not assume a Gaussian veloc- 
ity distribution and it represents the azimuthally averaged second- 
order velocity moment which includes rotation: 



2 1 ^ 



(2) 



where N is the sample size and Aui is the uncertainty on the radial 
velocity vt of the i* ^ globular clust er. The unc ertainty on ugc is 
estimated through the formula fiom lOanese et a l. ( 1980). We will 
refer to vgc of the red, blue and all GCs as vgc.,b., «gc,b, «gc 
respectively. 

The difference between ctgc and dgc is that the former repre- 
sents the rotation-subtracted velocity dispersion whereas the latter 
also includes the rotation of the spheroid and it is a better reflection 



of specific kinetic energy. A comparison between ctgc and ^gc is 
given in Figure[T]for our galaxy sample without any GC subpopu- 
lation split. The two quantities are consistent with each other when 
the rotation component is negligible, as seen for several systems. 

We perform a "sanity check" on all literature data. We prune 
GCs deviating more than ?)a from the local GC velocity distribu- 
tion. We also clip outliers with unreasonably large uncertainty (usu- 
ally > 100 km s^^) and then we recalculate the respective ctgc and 
VGC to avoid methodology biases. 



3.2 



The A/bh 
systems 



- ctgc and Afen — «gc relations for GC 



Here we describe how we characterize the A/bh — ctgc relation. 
The procedure is identical for the A/bh — t^GC relation. 

In logarithmic space, A/bh and ctgc appear to be linearly cor- 
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related. The relation we want to study is therefore: 



'°s(^)^" + ^^°§ 



0"GC 



200 km s" 



(3) 



where a and /3 are the intercept and the slope of the relation. 
The numerical constant (200 km s^^) is the normalization fac- 
tor adopted in similar studies of the stellar A fen — o-« relation . 
We then use t he y^-minimization te chnique l lPress et al.lll992h 
as modified by iTremaine et al.l ( l2002h . This ensures that the best 
fit to Equation [3] is not biased in the case of large uncertainties 
dParket al.ll2012h . Our minimization function is, using the nota- 
tion y = a + I3x: 



^K/?) = E 



{yi 



■Px, 



i4. + /3^<, + ^^ 



(4) 



where e^ and ey are the errors on x and y respectively. These 
are defined as e^, = (logfjuppcr - logcriower)/2 and ey = 
(log A/bh, upper - logA/BH,iowor)/2, respectively. The term e is 
the intrinsic scatter in the y direction in units of dex. e is iteratively 
adjusted so that the value of X^/i^ - 2) equals 1 ± %J2/N. Un- 
certainties on a and jS were obtained by bootstrapping the sample 
2000 times and selecting the 68 per cent confidence interval. 

This x^ estimator does not treat the data symmetrically in the 
presence of intrinsic scatter. An "inverse" regression (minimizing 
the log (TGC residuals rather than the log A^bh residuals) can lead 
to veiy different slopes. The latter is preferable in the p resence 
of possible Malmquist-type biases (see lGraham et alj201ll) . Given 
our ignorance of the physical mechanisms which links black hole 
mass to velocity dispersion, there is no reason to believe that the 
forward regression should be favored over the inverse regression. 
Therefore, we perform both the "forward" and the "inverse" re- 
gression by replaci ng e in Equation |4] with /3^e^ as suggested by 
lNovaketalJ ( l2006l) . 



4 RESULTS 

The A/bh — t^GC (and the A^bh — o"gc) diagrams for our sample 
are shown in Figure |2] in which the final slope and intercept of the 
relations are the average between the forward and the inverse fit. 
The respective best fit parameters are reported in Table [T] 

We find that A/bh correlates both with ace and ?;gc for all 
GC subsamples. However, we note that the intrinsic scatter of all 
our GC subsets are at least two times larger than those reported by 
SC12. We find that this disagreement is driven by the A'/bh val- 
ues of five galaxies in the SC12 sample (marked in Table lAl] with 
"a") for which we have updated A'/bh measurements. In fact, re- 
analyzing the SCI 2 sample using our new vel ocity dispersion val- 
ues and the A/bh values from SC12 (all from lGiiltekin et al .1120091 
and references therein), we always obtain e ^ 0.31 dex, which is 
in agreement with their findings. Conversely, the regression on the 
SCI2 sample using updated A/bh values, returns e = 0.38 dex for 
the full sample and e = 0.44 dex for the blue and red GC subsets. 
We conclude that the small intrinsic scatter of SC12 is driven by 
their black hole mass values and not by their GC system velocity 
dispersion data. This assumes that the latest values of A/bh that we 
adopt here are also more accurate than those which preceded them. 

The slope, intercept and intrinsic scatter of the A/bh — ugc 
and the A/bh — ctgc relations are consistent with each other within 
the errors. Similarly, the differences found for the blue and red 
GCs are not statistically significant. We note that the slopes of the 
A/bh — ^GC relations are always steeper than the A/bh — o"gc 
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Figure 3. Cumulative root-mean-squai'e velocity dispersion profiles. The 
plot shows the rotation-included velocity dispersion profiles within a certain 
radius for all (top panel), blue (central panel) and red GCs (bottom panel). 
A running mean is used. Different colors represent different galaxies. Most 
of the profiles are generally flat at all radii. 



ones, because iigc > o"GC at low masses. Also, the smaller intrin- 
sic scatter with vgc suggests that the GC kinetic energy (rotation 
plus dispersion) is a better predictor of black hole masses than the 
rotation-subtracted velocity dispersion. 

The intrinsic scatter of the A/bh — ''gc and A/bh — c"gc 
relations ar e slightly larger than that of t he stellar A/bh — o* re - 
lation from lMcConnell&MiJ ( 1201 2h and lCraham & Scotll ll2012h . 
who both find e ~ 0.4 dex. The best-fit to the stellar A/bh — cr« 
relation computed using our 21 galaxies has an intrinsic scatter of 
e = 0.35^0 Q5 dex in the log A/bh direction, which is also consis- 
tent with previous findings. 

Lastly, it is noted that the stellar A/bh — o"* relation in Figure 
|2]is shifted towards larger velocity dispersion values with respect to 
the A/bh— (GC system velocity dispersion) relations. This offset is 
expected because a, and the GC system velocity dispersion sample 
different regions, and maybe different physics, of the galaxy veloc- 
ity dispersion profile. The stellar velocity dispersion, which probes 
(/? < /?e), is usually larger than the GC system velocity disper- 
sion, which usually probes R > Re. The difference {a, — Wrms.A) 
is found to have a mean of 35 ± 6 km s^^ for our 21 galaxies. 



4.1 Radial trends 

It is interesting to see if the properties of the A/bh — I'GC or A/bh — 
(JGC relation vary when the velocity dispersion is computed within 
different galactocentric radii. 

To do so, we first normalize the galactocentric radii of each 
GC system to the host galaxy effective radius. We then perform 
X^ tests (Equation |4j with vqc and (jgc computed within dif- 
ferent radial bins. For the sake of consistency, we adopt effective 
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Figure 4. Best-fit Mbh — ''GC relation within different radial bins. The plot 
shows how the best-fit a, fi and e vary when the vqq is computed within an 
increasing number of effective radii. Blue and red colours represent the two 
GC subpopulations. SoHd and dashed lines are the results from the forward 
and the inverse regressions respectively. The hoiizontal axis is the radius of 
the outermost GC in a given radial bin. For clarity, only the error bars from 
the forward regression are shown. None of the radial trends are statistically 
significant. 



radii values from 2MA SS, and we use the transformations from 
ICappellari et al.l ( 1201 ll) to make them con sistent with the values of 
the RC3 catalogue l lde Vaucouleurslll99in . 

The cumulative velocity dispersion profiles for our galaxy 
sample are shown in Figure |3] for all GC subsets. The profiles are 
generally flat over the radial range probed. It is worth noting that 
GC dispersion profiles span different radial ranges depending on 
the galaxy, and we do not extrapolate the dispersion profiles to 
compensate for this effect. Therefore, the number of GC systems 
within a given effective radius varies with the radius itself. De- 
manding a minimum of six GC systems per radial bin, we study 
the Mbh — fGC and the A/bh — ctgc relations between 3.5 and 
5.5 Re for the blue and the red GC subpopulations. 

Results are shown in Figure |4] for the A/bh — fGC relation. 
Each radial bin contains between six to a maximum of eleven GC 
systems. As expected from the flatness of the velocity dispersion 
profiles (Figure |3), none of the radial trends seen in Figure 2] are 
statistically significant. The relations for the blue and the red GC 
subpopulations are also statistically indistinguishable. There is an 
hint that the intrinsic scatter for the A/bh — fGCR becomes smaller 
towards the central regions. This result is biased by the fact that 
the red GCs tend to be more centrally concentrated than the blue 
GGs. Given the small number statistics, the best fit to the A/bh — 
VGc relation is independent of radius within which the velocity 
dispersion is measured, at least for R > R^. The same exercise 
performed on the rotation-subtracted velocity dispersion ctgc leads 
to a similar result. 

A caveat to bear in mind is the way the GC system velocity 
dispersion is computed. Ideally, one should weight the velocity dis- 
persion for the GC surface density within a certain radius, similarly 
to what is done for the stellar v elocity dispersion cr* (see Equa- 
tion 1 in lMcConnell & Mal2012h . Similarly, the scale radius used 
in Figure [3] should be the GC system's effective radius and not the 
host galaxy's effective radius. However, GC surface density profiles 
are not available for all our galaxies. They are also dependent on 
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Figure 5. AfBH — ^GC relation. Orange and black points are galaxies 
with and without a core in the inner suiface brightness profile respectively. 
The filled and dashed lines are the best fits to core and non-core galaxies 
when using vqq and ctgc respectively. Green boxes mark barred galaxies 
(NGC 1316, NGC 1023 and NGC 253). The slope of the A^bh - '^GC re- 
lation for core galaxies is consistent within the errors with that of non-core 
galaxies. 



variables such as GC selection criteria and imaging field-of-view, 
which have been carried out differently in the literature. 

On the other hand, total GC system size scales with galaxy ef- 
fective radius (Kartha et al. in prep.) and we see no strong variation 
of GC system velocity dispersion with radius. 



4.2 Cores and bars 

The stell ar A/bh — o"* relation is different for galaxies with or with- 
out bars JGraham et al.ll20I If) . It is thought that the orbital struc- 
ture of the bar may elevate the apparent bulge velocity dispersion 
JBureau & Athanassoulalll999h . resulting in an offset A/bh — f7"» 
relation for barred galaxies with the appropriate bar orientation. On 
the other hand, the A/bh — cr* relation does not differ for non- 
barred galaxi es with or without a 'co re' in the inner surface bright- 
ness profile ( IGraham & Scottll20I2h . An exception may however 
exist for ult ramassive black holes s uch as those in NGC 3842 and 
NGC 4489 JMcConnell et al.ll201ll) . If these are included in the fit, 
the A/bh — cr. relation for core galaxies is steeper (/3 ~ 7.0) than 
that for non-core galaxies. 

We have tested if the trends seen for 'core' and barred galaxies 
with stellar data are also present in our A/bh — vac and A/bh — 
(TGC relations. To avoid low number statistics issues, we only look 
at the whole GC population, without any colour split. 

Our sample contains only three barred galaxies (NGC 1023, 
NGC 1316 and NGC 253), preventing any statistical analysis. For 
the sake of completeness, we note that NGC 1023 and NGC 1316 
are indeed offset to higher velocity dispersions relative to the best- 
fit A/bh — vac relation (Figure |5). However, only NGC 1023 is 
offset when considering ctgc- NGC 253 is neither offset from the 
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AfsH — fGC nor the J\/bh — o"gc relation, in agreement witli wliat 
was found for stellar data. 

Regarding 'core' galaxies, our sample contains nine core 
galaxies and twelve non-core galaxies (see Table IaH . The centre 
of the galaxy NGC 1407 is actually unclassified, but assume this 
galaxy to have a central core gi ven its mass. We treat NGC 1316 
(Fornax A) as a cored galaxy (IFaber et al.lll997h. but the re ader 
should see the cautionary remarks in lOraham & ScottI (120121) re- 
garding this galaxy's lack of a bulge/disc decomposition. 

The relation between A/bh and GC system velocity disper- 
sion for core/non-core galaxies is shown in Figure |5] We remind 
the reader that the final slope of the A/bh — (GC system velocity 
dispersion) relations is the average between the forward and the in- 
verse regression. Using the uncertainty on the slope and intercept 
of each regression, we derived a weighted mean to account for the 
large uncertainties caused by low number statistics. For non-core 
galaxies, we obtain a slope of /3 = 3.6 ± 1.5 and /? = 2.8 ± 1.5 for 
the A/bh — CTGC and Mbh — vqc relations respectively. For core 
galaxies, the uncertainty on the slope from the inverse regression is 
larger than the slope itself. This means that the final slope of this re- 
lation is driven only by that of the forward regression. In this case, 
we find ^ = 2.2 ± 1.6 and ,9 = 2.4 ± 1.6 for the A^/bh - ctgc and 
A/bh — «GC relations respectively. In conclusion, the A/bh— (GC 
system velocity dispersion) relations for co re and non-core galaxies 
are consistent with each other as found bv lGraham & Scotli ( l20I2r) 
with stellar velocity dispersion data. 
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[kms-1] 
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Table 2. Black hole mass predictions. Listed from the left to right ai'e: 
galaxy name, morphological type, GC root-mean-square velocity disper- 
sion, predicted black hole mass from Equation |5] and the GC refer- 
ences. Galaxies below the horizont al Une are from t he SLUGGS sur- 
vey. R efer ences to GC data are: 1. iPuzia et al.l J2004l); 2, iMisgeld et alj 
i201 ih: 3. iNorris et al.l 120121); 4. iLee et all bOlOh ri TlPark et alj J2012h: 
6. 7,lBeaslev et alj )2009hr8riBeaslev et alj )2006l) ; 9.10.11.12, |Pota et alj 
i2013h ; 13. lFosteretal.N20TT[) . 



4.3 Predicting AfBH in other galaxies 

We exploit the best fit A/bh— (GC system velocity dispersion) re- 
lations found in this work to predict A/bh in galaxies without di- 
rect black hole mass measurements. We collected a sample of 13 
galaxies with GC system kinematic inform ation, listed in Ta ble [2] 
The first four galaxies were re-analyzed in IPota et al.l (1201 3l) with 
the methods described in i ]3.1| S imilarly, we re-a nalyzed the GC 
system kinematics of NGC 4406 ( jParket al.ll2012l) and of thre e lu- 
minous Virgo dwarf ellipticals (dEs) from lSeaslev et al.l ( l2009l) and 
ISeaslev et alj |200d ). Given that the A/bh— (GC system velocity 
dispersion) relations for the blue and the red GC subpopulations 
return consistent results, we decided to use the best fit A/bh — wgc 
relation: 



log 



A/b 



Mq 



;.79 + 4.45 log 



«GC 



200 km s" 



(5) 



where the slope and the intercept of this relation are the average 
between the forward and the inverse regression from Table [T] 

Predicted black hole masses are given in Table |2] Particu- 
lar emphasis should be given to the three Virgo dEs, whose pre- 
dicted A/bh falls into the range of intermediate mass black holes 
(< IQp M(^). All thre e dEs are known to have a nuclear star cluster 
( JFerrarese et al.ll2006l) . whose masses are about one order of mag- 
nitude larger than our p redicted black hole masses, as is expected 
dScott & Grahamll20I3h . In fact, the relation between the mass of 
the nuclear star cluster A/nc and stellar velocity dispersion a, does 
not run parallel to the stella r A/bh — ct*. At fixed a, < 150 km 
s' MCraham & ScottI J2012h shows that A/nc > A/bh, which is in 
agreement with our findings. 

It is also worth noting that NGC 3311, the dominant ellip- 
tical galaxy of the Hydra Cluster, is at first glance predicted to 
host an ultramassive black hole candidate with A/bh ~ 10^^A/q. 
However, caveats here are the inclusion of ultra compact dwarfs 
(UCDs) which make up half of the kinematic sample of this galaxy. 



UCDs ca n be kinematically distinct from the underlying GC sys- 
tem (e.g. IStrader et al.ll20I lb and they can bias the velocity dis- 
persion calculation. Another source of contamination mig ht come 
from intra-cluster UCDs/GCs (Misgeld et al. 2011 ; Richtler et all 
uOIin . Excluding the 52 UCDs and looking only at the GC sam- 
ple, which may still be biased by the cluster potential, we obtain 
A/bh = 8.4125 x 10^^°. Thi s is still more massive t han the most 
massive SMBH known today iMcConnell et alj|201ll ). 



5 DISCUSSION AND CONCLUSIONS 

The aim of this paper was to test how well the velocity dispersion 
of extragalactic globular cluster systems correlates with the mass 
of supermassive black holes. This was motivated by the work of 
ISadoun & ColinI ( |20I2|) who found an intriguingly tight correlation 
using 12 globular cluster systems. 

In this work we have extended this study to a sample of 21 GC 
systems and we have used the latest compilation of SMBH masses. 
We confirm that the velocity dispersion of GC systems correlates 
with SMBH m ass. However, this corre lation is less significant than 
that inferred bv lSadoun & ColinI (|2012|) . The tight correlation found 
by these authors was driven by old, and possibly less accurate, 
black hole mass values. We observe an rms scatter in excess of 
0.4 dex in the log A/bh direction. 

We looked at the correlation between A/bh and the velocity 
dispersion of the blue and the red GC subpopulations separately. 
Different scatters are expected if blue and red GC systems trace 
the kinematics of the halo and the bulge of the host galaxy respec- 
tively. In the case of a stronger correlation with red GCs, this would 
suggest that SMBHs grew along with the stellar bulge. Conversely, 
a stronger correlation with blue GCs would suggest that SMBHs 
formation is more closely related with the halo. Our current results 
cannot discriminate between these two scenarios. In general, we 
find no significant difference between the A/bh — (GC system ve- 
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locity dispersion) relation for the blue and the red GCs. This can be 
due to some factors discussed below. 

Ideally, one should analyze the bluer and the redder GCs for 
each GC system to avoid contamination in proximity to the blue/red 
dividing colour. This can make a difference in th e final value of the 
GC system velocity dispersion jPotaet al.l2013h . At the same time, 
this would decrease the number statistics for most of the galaxies. 
Also, uneven GC spatial samplin g can affect the final k inematic 
outcome, as seen for NGC 4636 in lschuberth et al.l (l2012h . 

We have looked at the Mbh— (GC system velocity disper- 
sion) relation computing the GC system velocity dispersion within 
different galactocentric radii, obtaining no significant trends with 
radius. Collectively, this suggests either that the A/bh— (GC sys- 
tem velocity dispersion) relation is secondary, or that a larger 
galaxy sample will be needed to discriminate which of the GC sub- 
population trends is the stronger. 

We have looked for possible trends in the A/bh— (GC sys- 
tem velocity dispersion) relation for core/non-core galaxies, finding 
similar slopes, in agree ment with stellar velocity dispersion results 
JGraham & Scotlll2012h . 

The best fit relation between A/bh and the rotation-included 
GC system velocity dispersion has been used to predict black hole 
masses in 13 galaxies. This implies that NGC 3311 contains an 
ultramassive black hole with A^bh ~ 10^^ A/q. 
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14 
16 
+27 
25 
13 
13 



150 
152 
223 
162 
100 
148 
328; 
235 



17 
12 
15 
12 
9 
8 
9 
8 
+ 7 
6 

13 
10 
14 
12 
12 
11 

+ 12 

-9 



1451?? 
fi+21 



146: 



15 



23111? 
1661" 
1051^1 
13511? 
33711^ 
268l?l 



1541?? 
160l?6 

215llg 

1581- 

94l^« 

170l?^ 

2961^1 

203II2 



14 
18 
14 
15 
14 
16 
17 
14 
14 



Table Al. Galaxy samp le. Galaxy NGC names (1 ) and Hubble types (2) are from NED database. Galaxy d istances (3), SMBH masses (4) and stellar velocity 
dispersions (5) a re from lMcConnell & Mal l l2012[) and references t herein. If not in lMcConnell & Mal2012h : distances were obtained by subtracting 0.06 mag 
jMde^alJ 2007h to th e di s tance modulus from lTonrv et al.l 1 1200 ih : central stellar velocity dispersions are weighted values from HyperLeda; A/bh are from 
lOliva et al.l i 19951) and lHti 120081) for NGC 253 and NGC 5846 respectively. (6) is the presence of a core in the galaxy inner surface brightness profile. Column 



(7), (8) and (9) are the rotation-subtracted velocity dispersion for all, blue and red GCs res pectively. Column (10 ), (11) and (12) are t he root-mean square 
veloci ty f or all, blue and red GC s re spectively. GC references (13) : I.lLee et al.U2q08h : 2.l01senet alJJ2004|l : 3. lBeaslev et all j2qo4 : 4. iGoudfrooii et al J 
i200ll): S.ISchuberth et al.l <2010l); 6.lNantais & Huchral | |2010|); 7, jPierce 



11, 
17 
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<2010l); 6.lNantais & Huchral fcOlOl); 7, jPierce et al.ll2O06l): 8. JBergond et aUEooa); 9.|Puziaet alji2004l) ; lO. ICote et al.l i2003l) ; 
ll); 12. lHwang et al.l j2008l) ri3. IWoodlev et alj J2010l) : 14. |Pota et alj J2013h : I5. lAmold et alj J201ll) : 16. Foster et al. (in prep.); 



Strader et al.l J201ll); 18, Pota e t al. in prep. GC system velocity dispersion values for NGC 224 are from lLee et al.l 1 120081) . Galaxies with updated M] 
measurements after lGiiltekin et al.l J2009i) are marked with (a). 
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